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Using a generalization of the Poisson-Boltzmann equation, dynamics of counterion condensation 
is studied. For a single charged plate in the presence of counterions, it is shown that the approach 
to equilibrium is diffusive. In the far from equilibrium case of a moving charged plate, a dynamical 
counterion condensation transition occurs at a critical velocity. The complex dynamic behavior of 
the counterion cloud is shown to lead to a novel nonlinear force-velocity relation for the moving 
plate. 
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Many experimental techniques used to probe proper- 
ties of macroions involve setting them to motion by ap- 
plying some external driving mechanism (such as migra- 
tion in an applied electric field or hydrodynamic flow, 
sedimentation, etc.), and monitoring the dynamical re- 
sponse. A macroion in solution is surrounded by an ac- 
companying cloud of counterions, and the corresponding 
response of the counterions is in general very complex. 
This often makes interpretation of the observed prop- 
erties of the dynamically perturbed system a nontrivial 
task. 

For example, it is known that if a macroion moves, 
the asymmetry set by the motion causes a distortion in 
the counterion cloud. The distorted cloud, in turn, pro- 
duces a nonzero electric field (called the asymmetry field 
) that opposes the motion of the macroion. This phe- 
nomenon is known as the "relaxation effect," and has 
been studied in the literature within the Debye-Huckcl 
(DH) theory The approximate DH description, 

which is based on a linearization of the nonlinear Poisson- 
Boltzmann (PB) equation, is however known to break 
down when the macroions are highly charged. In partic- 
ular, it does not capture the phenomenon of counterion 
condensation, about which a lot has been learned by ex- 
amining the available exact solutions of the full nonlinear 
equilibrium PB equation |}J. 

In this Letter, we develop a dynamical generalization 
of the nonlinear PB equation, which we use to study the 
dynamical aspects of counterion condensation. For the 
specific example of a single charged plate in the pres- 
ence of counterions, we show that the equilibrium, cor- 
responding to a Gouy-Chapman profile, is approached 
diffusively In the far from equilibrium case of a moving 
charged plate we present an exact stationary solution of 
the dynamical PB equation, and show that the cloud of 
counterions phase separates into a comoving part and a 
condensate of the excess charge which is expelled to in- 
finity. The amount of the comoving charge depends on 
the velocity of the plate, and is decreased as the velocity 
increases. At some critical velocity v c , the system un- 



dergoes a dynamical phase transition, and eventually all 
the counterions evaporate from the neighborhood of the 
moving plate. We show that this leads to a novel nonlin- 
ear force-velocity relation for the moving charged plate 
in the presence of the counterions. 

To study the dynamics of counterions in the presence 
of macroions, one needs to formulate a dynamical gener- 
alization of the PB equation. Consider a set of negatively 
charged macroions described by the given charge density 
p m (x, t), embedded in a cloud of positively charged coun- 
terions with concentration c(x, t) 0] . The electric poten- 
tial ^(x,t) is then given by the solution of the Poisson 
equation: 
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V 4> = — (ec + p m ), 
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(1) 



in which e is the dielectric constant of the solvent, and e is 
the electron charge. The dynamics of c(x, t) is governed 
by a continuity equation of the form dtc + V • J = 0, 
in which the current density J is composed of a deter- 
ministic part cv, and a stochastic part — DVc as given 
by Fick's law, where v(x, t) is the velocity field, and D 
is the diffusion constant of the counterions. In a mean- 
field approximation, the velocity field at each point is 
determined by the local value of the electric field E(x, t) 
as v = /iE = —fiV(f>, where /i is the electric mobility 
of the counterions. Using the self-consistency relation in 
the continuity equation, one obtains the so-called Nernst- 
Planck equation 



d t c 



DV 2 c + ^V ■ (cV0). 



(2) 



Note that we have simplified the problem by neglect- 
ing any couplings to the hydrodynamics of the solvent. 
The above two equations (Eqs.(|l|) and (^)) describe the 
nonlinear dynamics of counterions in the presence of a 
given distribution of macroions. At equilibrium, d t c = 
and Eq.(^) yields c ~ exp(— fj,<p/D), where equilibra- 
tion at a temperature T implies an Einstein relation 
D = (fi/e)kBT. Inserting the Boltzmann form for c back 
into Eq.(ffl) then yields the PB equation. In analogy to 



1 



the equilibrium case, we can define a dynamical Bjerrum 
length Ib = ep,/eD. 

The above nonlinear equations are in general very dif- 
ficult to solve, which is why they are usually dealt with in 
a linearized approximation |^|,^). However, the approx- 
imation is known to break down in the neighborhood 
of highly charged surfaces, where interesting phenomena 

d t <j>(x,t) - d t 4>{-oo,t) = Ddlcj>+^{d x <l>f + — f 



such as counterion condensation take place. To capture 
the essence of the nonlinearity in the dynamical context, 
we restrict ourselves to a one dimensional case where the 
equations prove to be more tractable. Eliminating the 
concentration field c from the two equations, one obtains 
an equation for the electrostatic potential, which can be 
integrated twice to yield 

[^p m (x',t)d x (j)(x',t) + Dd x p m (x',t) + J m {x',t)] , (3) 



where the macroion current density J m is obtained from 
the continuity equation for the macroions dtp m +d x J m — 
0, and the boundary condition E(—oo, t) = has been 
implemented. 

The above equation, which can be called a dynamical 
PB equation 0], belongs to the general class of the cel- 
ebrated Kardar-Parisi-Zhang (KPZ) equations [||, with 
the coupling to the source appearing both in the additive 
and multiplicative forms. It is however important to note 
that in this context the sources represent slower dynam- 
ical degrees of freedom corresponding to the motion of 
the macroions (compared to counterions), which makes 
it somewhat different from the standard KPZ problems 

To study how counterions dynamically rearrange them- 
selves around macroions, we focus on the specific ex- 
ample of a single negatively charged plate with charge 
(number) density a and surface area A, which is moving 
at a constant velocity v, i.e., p m {x,i) = —e<j5(x — vt). 
Using a Cole-Hopf transformation (j)(x,t) — 4>(— oo,t) = 



2D/ f i \nW{x,t) §,0], Eq.(|) can be written as 
d t W = DdlW 



2D 



5{x — vt) 



v - pE 
D 



e(x- vt) 



W, (4) 



where Eo(t) = —2D/p, d x In W\ x=v t is the asymmetry 
field, &(x) is the step function, and A = 1/ttoIb is a 
dynamical Gouy-Chapman length. Note that despite its 
simple form, this diffusion-like equation is still nonlinear 
due to the self-consistent coupling of the asymmetry field 
(electric field at the charged plate). 

Let us first assume that the plate is not moving (v = 0), 
and ask how the equilibrium configuration is dynamically 
approached. If we assume that the initial configuration of 
the counterions is symmetric with respect to the plate, 
the symmetry will be conserved during time evolution 
rendering Eo(t) = at all times. This simplifies Eq.(Q) 
to a linear diffusion equation which can be solved exactly. 
One obtains 
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In simple terms this means that the equilibrium Gouy- 
Chapman profile diffusively develops, as the initial con- 
figuration diffusively disappears. Naively, one might not 
expect such an amazingly simple behavior from a fully 
interacting Coulomb system. 

For a stationary plate at equilibrium, all the counte- 
rions are known to be "condensed" in a Gouy-Chapman 
(density) profile which decays algebraically at large sep- 
arations (see below) . It is interesting to study this prob- 
lem in the far from equilibrium case of a moving plate. In 
particular, we may ask if there is a stationary solution of 
Eq.(||), of the form W(x,t) = W(x — vt), in which a^the 
counterions are comoving with the plate at the velocity 
v. The answer to this question turns out to be negative. 
Apparently, the Coulomb attraction which could over- 
come the combination of entropy and repulsion among 
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the counterions in the equilibrium case, is not capable 
of accommodating the viscous drag of all the counteri- 
ons in this balance, for any nonzero v. However, we can 
find stationary solutions if we do not require that all the 
counterions are comoving. 

Such solutions can be shown to exist only for neg- 
ative values of v, with a boundary condition that re- 
quires a specific value for the electric field at +oo. This 
implies that the neutralizing counterion cloud actually 
phase separates into a comoving domain and a domain 
of excess charges accumulated at +oo, as required by 
this one dimensional model This class of solutions 
is parametrized by the number of comoving counterions. 
The stationary solution of Eq. (Q) with maximum number 
of comoving counterions can be obtained as |9] 



1 - 



W(x,t) = 



1+A"/4D \ 
l-\v/AD ) eX P 



v(x— vt) 
D 



(l-Xv/AD)' 2 \ 



1 



l-(^) 2 ]^}ex P 



v{x—vi) 
215 



for x < vt, 



for x > vt, 



(6) 
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and correspondingly for the counterion density profile as 



c(x,t) = < 



2tt£ b 



1 • 

\, /l + A^/4D\ 2 f v(x-vt)l I 



1 \l-(\v/±D) 2 } 



for x < vt, 



for x > vt. 



(7) 



The above solution has very interesting features. In 
the limit v — ► 0, Eqjffl) reduces to the equilibrium Gouy- 
Chapman profile ]10| 
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In this limit, all the counterions are condensed, typically 
at a distance given by A, although the profile has an al- 
gebraic decay. 

For a nonzero v, while the density profile is still alge- 
braically decaying "behind" the moving plate (for x > 
vt) , it decays exponentially "ahead" of the moving plate 
(for x < vt) , indicating a jamming of the profile as a re- 
sult of the motion. In the comoving profile of counterions, 
there is a total charge of Q> = |Qo[l — (Xv/AD) 2 ] be- 
hind the plate, as opposed to Q< = |Qo[l + (Xv/AD) 2 + 
Xv/2D] ahead of the plate, where Qo = eaA is the overall 
charge of the counterions. Note that the comoving coun- 
terions are distributed asymmetrically with Q> > Q<. 
The amount of the excess charge can be easily obtained 
as 



Qeyi(v) = Qo - Q> - Q< 



eA 



D ' 



(9) 



which is independent of ex. By examining the asymptotic 
limit of the electric field at infinity, one can easily see 
that this excess charge is indeed expelled to +oo. 

As the velocity of the plate increases, the number of co- 
moving counterions decreases, until at Xv/AD — — 1 even- 
tually all the counterions evaporate to infinity. This cor- 
responds to a dynamical counterion condensation phase 
transition, happening at the critical velocity v c — 4D/A, 
and is similar to the equilibrium transition in the case of 
a charged cylinder pl| , ^2| . It is easy to check from Eq. (Q) 
that the density profile of the comoving counterions van- 
ishes on both sides at Xv/AD = — 1. 

The solution of Eq.(^) yields a value Eq = (1 + 
Aw/8-D)u//x for the asymmetry field at the position of 
the plate, which leads to a very interesting mechanical 
response for the charged plate |Tj| . In the stationary sit- 
uation, the total force F to t exerted on the charged plate 
should vanish. This means that an externally applied 
(mechanical) force F ext (that is necessary to maintain 
the constant velocity motion) should balance an electri- 
cal contribution due to nonzero value of the electric field 



Eq felt by the charged plate, namely, F cxt — QqEq = 0. 
Inserting the velocity dependent value for Eq, one obtains 



l-F«xtl = 



(Qo/aO \v\ - ^ (v/fJ,) A, for \v\ < v c , 



2irQ 2 /eA, 



for \v\ > v c , 



(10) 



as the force-velocity relation for a moving charged plate. 
A plot of Eq.([l0|) is sketched in Fig. 1, in which F c — 
2irQ 2 /eA. 
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FIG. 1. A plot of nonlinear force-velocity curve for 
the moving charged plate. The slope of the curve at any 
velocity is a measure of the number of condensed (co- 
moving) counterions. 

In the limit of very low velocities (\v\ <C v c ), one can 
neglect the quadratic velocity term in Eq. (|To| ) . This leads 
to a viscous drag with an effective friction coefficient 
Ccff = Qo/fJ-, where the "electrostatic friction" is due to 
dissipation corresponding to the comoving could of coun- 
terions. Note that all the counterions are condensed in 
this limit, and participate in the dissipation. In the op- 
posite limit of \v\ 3> v Cl no counterions are condensed 
and the friction coefficient is consequently zero. How- 
ever, since all the counterions participate in forming a 
condensate at infinity, there is a constant (independent 
of velocity) electrostatic attractive force acting on the 
plate, which is simply the force between two parallel and 
equally charged plates. In fact, one can rewrite Eq.dTTJ) 
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in such a way that these features are manifest at any 
velocity: \F ext \ = (Q c (v)/ fi)\v\ + 2nQ c ^(v) 2 / eA, where 
Q c (v) = Q > + Q < — Qo — Qcx{v) is the amount of con- 
densed (comoving) charge. Interestingly, the nonlinear 
response function 



C{v) = 



_ dF ext (v) 



dv 



Q c (v)/n, 



(11) 



gives a proper account of the number of condensed (co- 
moving) counterions at any velocity. 

The dynamical counterion condensation phase transi- 
tion discussed above can be characterized by an order 
parameter, which is given by the number of condensed 
counterions. The tunning parameter for this dynamical 
transition is the external force, which is the analogue of 
temperature in the equilibrium case. In analogy to equi- 
librium critical phenomena, one can look at the critical 
exponents corresponding to the transition. For example, 
the order parameter vanishes at F c as Q c ~ (F c — |i 7, ext |)' 3 , 
with a mean-field exponent /3 = 1/2. 

It is well known that at equilibrium counterion conden- 
sation is determined by dimensional considerations p[ . 
In d dimensions, a D dimensional macroion of size Lu at- 
tracts a counterion at a perpendicular distance L± , with 
a Coulomb energy that goes like E c ~ 1/ L d ^ D ~ 2 . On 
the other hand, entropy of a particle confined in such 
a box can be obtained as S ~ ln(£j? 'L d ^ D ). Compar- 
ing energy and entropy one obtains that for d < D + 2 
counterions tend to condense to minimize energy (lower 
values for L± are favored), while for d > D + 2 they 
prefer to be free and gain entropy (higher values for L± 
are favored). The case d = D + 2 is marginal, where 
a counterion condensation transition takes place plff2| . 
The analysis presented in this Letter clearly shows that 
this simple equilibrium picture is modified when the sys- 
tem is far from equilibrium. In particular, we found that 
complete condensation in the case of a charged plate will 
no longer hold for a moving plate. Another important 
difference is that, unlike the equilibrium case, there is 
a marked contrast between the condensed and the free 
counterions. 

One might question the validity of the present formula- 
tion at high velocities such that Xv/D ~ 1, because of the 
assumption that the dynamics of the macroions is much 
slower compared to the counterions. The dimensionless 
parameter Xv/D (often called the Peclet number) can be 
written as r c /r TO , in which r c is the time it takes a coun- 
terion to diffuse a distance A, whereas r m is the time in 
which the macroion drifts a same distance. It is thus pos- 
sible that this ratio becomes comparable to 1, while the 
ratio between the corresponding diffusion times is still 
much less than unity. 

In the derivation of Eq.(ph, we have made use of a 
mean-field approximation similar to what is used to ob- 
tain the PB equation. It is well known that PB equation 
is incapable of describing intriguing phenomena such as 



attraction between like charged objects in the presence 
of condensed counterions, because of this approximation 
]i~5| |. One should then attempt to go beyond mean-field 
theory and study the effect of fluctuations on the dynam- 
ical PB equation along the lines of what has been done 
in the equilibrium case Jl6|-^8[|. In particular, it is shown 
in Ref. |Di] that fluctuations lead to a stronger condensa- 
tion, because mean-field theory overestimates the repul- 
sion between counterions. It will thus be interesting to 
see how this effect competes with the viscous drag in the 
case of a moving plate. Another important and very in- 
teresting extension of this work would be to incorporate 
the coupling to hydrodynamics of the solvent, which we 
anticipate to have dramatic effects. 
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